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^y^ \ The use of leading order effective field theory (EFT) to describe neutron- 

iy-\ ' deuteron scattering leads to integral equations that have unusual behaviour: 
when only two-body interactions are included, the scattering amplitude does 

'►-j^ . not approach a limit when the cutoff used to solve the equations is removed. 

r — ! It has recently been shown that this cutoff dependence can be eliminated by 

^^ \ the careful inclusion of a three-body force. In this paper we show that the 

Q^ ' cutoff dependence is just a reflection of the fact that the aforementioned inte- 

o 

f^ . corresponds to the physical scattering amplitude. We show how to numeri- 

^^ \ cally extract the physical scattering amplitude from the general solution and 

in this way explicitly demonstrate that the amplitude for a particle scatter- 



^ 

^ 



X 



Q ' ing off a two-body bound state, in leading order EFT, is in fact determined 



entirely by two-body forces. 
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I. INTRODUCTION 

Effective field theory (EFT) approaches to problems in nuclear physics have been under 
intensive investigation during the last few years. A review of recent developments (and 
references to the relevant papers) can be found in |jl|. In the two-body sector, most of the 
theoretical problems encountered initially have now been resolved with recent EFT calcula- 
tions of nucleon-nucleon scattering giving a particularly successful description [0-^. On the 
other hand, application of the EFT program to the three-body sector remains problematic. 
Here we would like to address an especially intriguing problem encountered in the lowest 
order EFT calculation of neutron-deuteron scattering in the J = 1/2 channel 0. 

For very low energy EFT calculations of nuclear systems, one integrates out all parti- 
cles other than the nucleon, and at leading order one is left with just a constant contact 
interaction between the nucleons. A problem arises immediately when one tries to use this 
contact interaction in a standard calculation of three-body scattering. In this regard we 
note that the spinor nature of the nucleons plays no essential role here and that exactly the 
same problem occurs in the case where the three particles are all scalars. Thus for the sake 
of simplicity and without loosing generality, we shall present our discussion for the case of 
three bosons, described by EFT, where at very low energies all other particles have been 
integrated out and the leading order interaction is a constant contact term. 

The problem in question can be described as follows. To find the leading order amplitude 
for a boson scattering off a two-boson bound state we need to sum up an infinite number 
of diagrams involving the two-body constant contact interaction. While each leading order 
three-body diagram with re-summed two-body interactions is individually finite, when one 
tries to sum up these diagrams by using an integral equation 

/■oo 

a{p,k) = M{p,k)+ dq Z{p,q)a{q,k), (1) 

Jo 

where a is the summed amplitude, M is the Born term, and Z is the kernel, one finds 
that the operator (1 — Z)~^ does not exist. Thus one cannot express the amplitude as 
a = (1 — Z)^^M with the consequence that the numerical solution of Eq. (|I]) is particularly 
difficult to construct. 

If one tries to handle this problem by introducing an ultraviolet cutoff for the integral in 
Eq. (|ID, one then finds that although the modified equation is easily solved (e.g. by matrix 
inversion), its solution is sensitive to the chosen cutoff and the limit as the cutoff is removed 
does not exist. To resolve this problem, Bedaque et al. introduced a one-parameter three- 
body force counter-term into their leading order EFT calculations. Indeed, they argued that 
the introduction of this three-body force is a necessary and sufficient condition to eliminate 
the cut-off dependence. 

In the present paper we solve the problem of sensitivity to the ultraviolet cut-off without 
introducing three-body forces into the leading order Lagrangian. We do this by recognising 
that Eq. (|T]) has an infinite number of solutions, only one of which corresponds to the 
physical amplitude for particle-bound state scattering. Furthermore, we show that the 
introduction of a particular ultraviolet cutoff into Eq. (|1|) results in an equation whose 
solution is an approximation to just one of the infinite number of solutions of the original 
Eq. (ffl). Moreover, by varying the ultraviolet cutoff, we end up obtaining approximations 



that jump between different solutions of Eq. (p. It is this jumping between solutions which 
is responsible for the sensitivity to the ultraviolet cut-off observed in Ref. 0. By performing 
a straightforward numerical analysis of the solutions with different cutoff parameters, we are 
able to construct the actual physical solution of Eq. (JlD where the input consists of two-body 
interactions only. In this way the present paper extends the ideas presented in Ref. 0] and 
also complements these through explicit numerical calculations. 

Although we have considered the case of three boson scattering for simplicity, exactly the 
same considerations hold for the problem of J = 1/2 channel neutron-deuteron scattering 
in the framework of EFT. We will address the specific case of neutron-deuteron scattering 
in a separate paper. 

II. THREE-BODY SYSTEM IN LEADING ORDER EFT 

A. Integral equations for boson-dimeron scattering 

The leading order Lagrangian for the considered EFT of non-relativistic self-interacting 
bosons is given by 

/: = 0^(^5o + ^)0-^(0W (2) 

where is the boson field, m is its mass, and Co is a coupling constant. For the sake 



of convenience |T0] one can rewrite this theory by introducing a dummy field $ with the 



quantum numbers of two bosons, referred to as a "dimeron" [^: 

£ = 0t Lq^ + Z_ ] + A$^$ - 4=('^V0 + h.c). (3) 

V 2m/ V2 

The scale parameter A is included to give the field $ the usual mass dimension of a heavy 
field. Observables depend on the parameters of Eq. (j^) only through the combination Co = 

The (bare) dimeron propagator is a constant i/A and the boson propagator is given 
by the usual non-relativistic expression i/{p^ — p^ /2m + it) where p = |p| (similar three- 
vector notation is used below for other momentum variables). The dressing of the dimeron 
propagator is illustrated in Fig. |l](a). Summing loop-diagrams, subtracting the divergent 
integral at P" = P = (where P° and P = |P| refer to the momentum of the dimeron) and 
removing the cut-off, one gets the following dressed dimeron propagator 0: 

^S{P\ P) = ^p (4) 

_A« + ^J-mP° + PV4 -te + ie 

where A"^ is the renormalised parameter (A has absorbed the linear divergence). Attaching 
four boson lines to this dressed dimeron propagator one gets the two-particle scattering 
amplitude at leading order. This amplitude has the form of an effective range expansion 
truncated at leading order with g^/A^ = ATia2/m where 02 is the two-body scattering length. 
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(b) 
FIG. 1. (a) Dressing of the dimeron. (b) Diagrams contributing to particle - dimeron scattering. 



For the scattering of a particle off a two-body bound state, standard power counting 
shows that the leading order contribution to the amplitude T is given by the diagrams illus- 
trated by the first equality of Fig. |l|(b). The sum of all these diagrams satisfies the equation 
represented by the second equality in Fig. |l|(b). For s-wave scattering it is convenient to 
define the function a{p, k) in terms of the s-wave amplitude To{p, k) by 



a(p, k) 



p^-k^ 



To{p,k) 



which can be shown to satisfy the equation |Tl|-[l4| 



(5) 



a{p, k) = M{p, A;) H / dq M{p, q) 

-^ Jo 



q' 



q^ — k? — ie 



a{q,k), 



(6) 



where k (p) is the incoming (outgoing) momentum magnitude, E = Sk'^/Am — l/ma\ is the 
total energy, and 



6 \a2 V 4 I pq \q^ — qp + p'^ — mh 



(7) 



Eq. (^) was first derived by Skorniakov and Ter-Martirosian |Tl| (S-TM equation) and has 
A = 1 for the three-boson case. Three nucleons in the spin J = 1/2 channel obey a pair of 
integral equations with similar properties to this bosonic equation, and the spin J = 3/2 
channel corresponds to A = —1/2. 

It was shown by Danilov |T^ that for A = 1 the homogeneous equation corresponding to 
Eq. (D, 



2A f°° Q^ 

ah{p,k) = — dqM{p,q)— -ah{q,k), 

TT Jo q'^ — k"^ — te 



(8) 



has a solution for arbitrary E. In particular, there exists a solution of the homogeneous equa- 
tion for every energy corresponding to the scattering of a projectile off a two-body bound 
state. Although all such scattering energy solutions are unphysical, it must be emphasised 
that they are purely an artifact of having two-body effective potentials that are zero-range 
((5-function potentials in coordinate space). By contrast, the solution of the correspond- 
ing homogeneous equation for non-zero range potentials (namely the bound state Faddeev 



equation) has no solutions for energies corresponding to scattering. In turn, it should be 
remembered that zero range effective potentials are themselves an artifact of restricting the 
EFT model to the lowest order terms - they are not a property of the full EFT approach [Q. 
Although unphysical, the existence of a scattering energy solution ah of Eq. (||) has 
practical consequences for finding the physical scattering amplitude aph satisfying Eq. (^. 
The problem is that the existence of ah implies that Eq. (^ has an infinite number of 
solutions given by a = aph + Cah where C is an arbitrary parameter. As we shall see, Eq. (^ 
has actually more than one solution for any given E. Writing these solutions as a^^ where 
i = 1,2,3,..., the most general solution of Eq. (^) can therefore be written as 

a = aph + ^ Cia\ . (9) 

i 

Thus the sum of the diagrams in Fig. |I](b), which defines the physical amplitude aph, is only 
one of an infinite number of solutions to Eq. (^. 

Setting A = 1 and writing Eq. (^) and Eq. (^) in operator form as 

a = M + MGoa (10) 

and 

ah = MGoah (11) 

respectively, it is clear from the existence of a non-zero solution ah of Eq. (|TT]) that the 
inverse operator (1 — MGo)~^ does not exist. This in turn means that one cannot write 
the solution of Eq. (|l^) as a = (1 — MGq)'^M with the practical consequence that Eq. (|[) 
cannot be solved by matrix inversion. Indeed we find that many other numerical methods, 
for example Fade approximants, are likewise unstable for this case. 

The task of finding the physical amplitude aph appears to be formidable. Even if Eq. (^ 
could be solved, one would still need to determine the appropriate values of the parameters 
Gi in order to extract aph- Attempts to find aph by writing it directly as the sum 

aph = M + MGqM + MGoMGqM + ... (12) 

also do not help, as this sum cannot be evaluated numerically due to extreme sensitivity to 
roundoff errors. Again, this numerical instability appears to be linked to the non-existence 
of (1 — MGo)~^. Fortunately, all these numerical difficulties can be overcome. The rest of 
this paper is devoted to accomplishing this task. 

B. Extracting the physical boson-dimeron amphtude 

1. Asymptotic behaviour 

In order to distinguish the physical boson-dimeron amplitude aph from the infinite number 
of non-physical solutions given by Eq. @, it is useful to examine the asymptotic behaviour 
of the general solution a{p, k) to the S-TM equation for large p. It has the form |TB 



a{p,k) = J2Adk)p'' + 0il/p) (13) 



where Si are roots of the equation 

8A sin7rs/6 



v^scos7rs/2 



0. (14) 



The summation in Eq. ( pTSD goes over all solutions of Eq. ( p!4D for which |Res| < 1. For 
A = 1 Eq. (p!^ has two roots for which |Res| < 1: s = ±iso, where Sq ~ 1.00624, so that 
Eq. (|13D gives the asymptotic behaviour of the boson-dimeron amplitude as 

aip,k)r^A,{k) p''^ + A2{k) p-'^o . (15) 



It can be shown [|T^ that for A = 1 the asymptotic behaviour of the solution a/i(p, k) of the 
homogeneous equation, Eq. (|), is also given by the right hand side of Eq. (|15|). 

What distinguishes the physical amplitude aph from any other amplitude a satisfying the 
S-TM equation is that it is the sum of the diagrams in Fig. |](b); that is, it is equal to the 
sum of the series obtained by iterating Eq. (H). In this way aph can be expressed as a power 
series in the parameter A. By contrast, it can easily be seen that for the physical amplitude 
Opft, the right hand side of Eq. (pTsP must be identically zero, i.e., it cannot be expressed as 
a non-trivial power series in A.0 Thus 

O'phiPi k) -^ Q as p — i> oo. (16) 

2. Unitarity 

Another essential property of the physical amplitude aph is that it satisfies unitarity; that 
is, the on-shell physical amplitude a = aph satisfies the relation 

a{k, k) - a*{k, k) = 2iXk \a{k, k)f . (17) 

In this respect it is important to note that if an amplitude a satisfies the S-TM equation, it 
does not necessarily mean that a satisfies unitarity. Indeed, if we write the S-TM equation 
in operator form as a = M + MG^a and then try to prove unitarity in the usual way, we 
would firstly want to write M~^ = a^^ + Gq, then subtract the Hermitian conjugate of this 
equation, and lastly, rearrange to obtain the unitarity relation a — a'^ = a{Go — GQ)a'^ which 
reduces to Eq. ( p!7D on-shell. However, if the inverse [I — MGq)^^ doesn't exist, then it is easy 
to show that either operator a has no right-inverse or M has not left-inverse. In either case, 
the usual proof of unitarity breaks down at the first step. On the other hand, because aph can 
be expressed as an iteration of the S-TM equation, aph = M + MGqM + MGqMGqM -|- . . . , 
the unitarity relation for aph can easily be proved directly without the need to take the 
inverse of any operator. 



^For if Ai and A2 were non-vanishing and we could write so(A) = cq + ciA + C2A^ + . . . , then 
substituting this into Eq. ( [TsD we would obtain a power series in A whose coefficients do not have the 
same asymptotic p behaviour as the corresponding diagrams coming from the iteration of Eq. (P). 
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As the physical amphtude satisfies unitarity, we would like to consider only those solu- 
tions a of the S-TM equations that likewise satisfy unitarity. For Eq. (^, this means a strong 
restriction on the allowed values for the parameters C,. This restriction on the solutions of 
Eq. (^ can be achieved through the introduction of the K matrix equation 

2A 1"°° a^ 

K{p, k) = M{p, k) + —P dq M{p, q)^^K{q, k), (18) 

TT Jq q — k 

where P stands for "principal value" and where only the real solutions for the K matrix 
K{p, k) are considered. It is then easy to show that the amplitude 



satisfies both the S-TM equation and the unitarity relation of Eq. (|r?|) . Conversely, any 
amplitude a(p, k) that satisfies the S-TM equation and is unitary can be shown to be of the 
form given by Eq. ([I9|) where K is real and satisfies Eq. ([T8|) . 

One might hope that the unitarity condition would single out the physical solution from 
an infinite number of solutions of Eq. (P). This would require Eq. (|I8|) to have a unique 
solution corresponding to the physical amplitude. Unfortunately the homogeneous equation 
corresponding to Eq. ([T8|) admits non-trivial solutions, and this means that Eq. (0) has in 
fact an infinite number of solutions. Nevertheless, we know from the above discussion of the 
asymptotic behaviour of a^^ and the relation between the a amplitude and the K matrix, 
Eq. ([T9|) , that amongst all the solutions to Eq. (plSf), the one that corresponds to the physical 
amplitude is the solution which vanishes for large p. 

III. NUMERICAL SOLUTION 

A. Solution to the inhomogeneous equation 

As we are specifically interested in the description of three bosons within leading order 



EFT, we shall implicitly assume that A = 1 in all the equations below. Writing Eq. ( p!8[ ) 
symbolically for this case as A' = M + MGqK, the fact that the homogeneous equation 
corresponding to this equation has a solution for an arbitrary value of the energy means 
that the inverse operator {1 — MGq)~^ does not exist. Thus trying to solve Eq. (p!8D directly 
presents exactly the same type of numerical difficulties as discussed above for the case of 
Eq. (^. In order to solve Eq. (|18D numerically we introduce regularisation such that, in 
contrast to Eq. (|18|), the regularised equation has a unique solution. Although this solution 
depends strongly on the cutoff parameter(s) of the regularisation, we will see that this is not 
an indication of any difficulty with leading order EFT, as suggested in Ref. [0], but rather 
it is simply a consequence of the fact that different regularisations correspond to different 
solutions of the unregularised equation. 
We consider the regularised equation 

2A f^ a"^ 

K\p, k) = M{p, k) + —P dq M\p, q)^^K\q, k) (20) 

IT Jo q^ - k^ 
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In (pa^) 
FIG. 2. The dot-dashed, dashed, and dotted hnes are solutions K '{p,k) of the inhomoge- 



neous K matrix equation, Eq. (pO]), corresponding to three different cutoffs Ai = 1.6 x 10^ i 



-1 



A2 = 3.2 X lO^a^ and A3 = 6.4 x lO'^a^ , respectively. The solid line is the function 
Kh{P) k) = K^^{p, k) — K^{p, k) which forms a solution of the corresponding homogeneous equa- 
tion. In all cases k = 0.82a7 . 



where 



M^(p,g) = M(p,g) + 



A + g2 



(21) 



M^G^] 



^ exists and Eq. (^) is easily 
K^, satisfies Eq. ( p!8|) with very 



For this regular isat ion the inverse operator (1 
solved. We find that the solution of this equation, K 
good accuracy for any A that is sufficiently large; moreover, by increasing A one can ob- 
tain solutions K'^ that will satisfy Eq. ([T8|) to any given accuracy. As expected from the 
asymptotic behaviour specified by Eq. (|T3|), we find that K^{p, k) has oscillating behaviour 
for large p - see Fig. ^ We also find that K^{p, k) has an oscillating behaviour with respect 
to A, in agreement with what was observed in Ref. [^. 

More information on the functional form of K^{p,k) can be obtained by constructing 
solutions K^*{p, k) (? = !,... ,4) corresponding to four different values of A. As each of the 
K^^{p, k), for Aj large enough, is a solution to Eq. (18), the differences K^*{p, k) — K^^{p, k) 
are clearly solutions to the homogeneous equation corresponding to Eq. (|l^). However, what 
is particularly interesting about such differences is that their ratio 

K^^{p,k)-K^^{p,k) 
K^^{p,k)-K^i{p,k) 

is found to be totally independent of the momentum variable p. From this observation we 
conclude that the structure of K^ must be of the form 



K^'ip, k) = Ko{p, k) + C{K)Kh{p, k) 



(22) 



where Ko{p, k) is a solution to Eq. (|TB|), Kh{p, k) is a solution to the corresponding homoge- 
neous equation, and C(A) is purely a function of A (we consider the initial momentum k to 



be fixed). As the normalisation of tlie function Kh{p, k) is not determined by the homoge- 
neous equation, for definiteness we take Kh{p, k) = K^^{p, q) — K^'^{p, q) which is drawn in 
Fig. together with three of the particular solutions to the inhomogeneous equation. Note 
that KQ{p,k) and Kh{p,k) do not depend on A. Also, because each K^{p,k) oscillates as 
a function of p with a phase that depends on A, Kq{p, k) and Kh{p, k) must similarly oscil- 
late but with differing phases. Being a solution to the inhomogeneous K matrix equation, 
Ko{p, k) can be written generally as Kph{p, k) + Kl{p, k) where Kph is the physical K matrix 
which has the vanishing asymptotical behaviour (i.e. the K matrix corresponding to the 
physical amplitude aph), and Kl{p, k) is another solution to the homogeneous equation. We 
may thus write the solution to Eq. (|1^) for large A as i^ = K^ where 

K\p, k) = Kpf^ip, k) + Kl{p, k) + C(A)ir,(p, k). (23) 



For large p the structure of Eq. ( pS] ) agrees with the asymptotical behaviour of the solution 
to the regularised equation obtained in 0:0 

K'^ip » l/a2, k) = ^-^ —— cos f So In ^\ (24) 

cos [sq In [p^a2) + ^ \ P* J 

where p* = exp {—5/sq) A and 5, 7 and e are cutoff-independent constants. Indeed, by 
writing Eq. (p3) as 



K^{p ^ 1/(325 k) = —7 cos [sq In (^02) + e] — 7 tan [sq In (Aa2) — (5 + e] sin [sq In (^02) + e] 

(25) 

and noting that the physical K matrix Kp^ vanishes for large p, we can deduce 
that Kl{p,k) ~ cos [soln (pa2) -l- e], Kh{p,k) ~ sin [sqIu (^02) + e] and C(A) ~ 
tan [soln (Aa2) — 6 + e]. The deduced asymptotic behaviour of Kh and K^ is consistent 
with Eq. (|l^) and is further borne out by our numerical results. 

The existence of (at least) two linearly independent solutions to the homogeneous equa- 
tion may seem a little surprising in light of the fact that Eq. ( [2^ ) contains just one solution, 
Kh, with a A-dependent coefficient. To further check that the structure given by Eq. (^) is 
consistent with there being more than one linearly independent solution to the homogeneous 
equation, we consider the equation 

2A r°° 

k{p, k) = M{p, k) + — P dq M{p, q)k{q, k), (26) 

71" Jo 



where 



M(p,i) = AiiJ4±Pi±P!|. (27) 

V3q W-pq + P^ ' 



^Note that the cutoff used in Ref. Q is different from ours; however, the functional form of the 
asymptotical behaviour does not depend on the particular choice of the cutoff. 



To obtain M we substituted a2 = oo and niE = in the expression for M{p,q), and 
also removed q^ /{q^ — k"^) from the integrand. The homogeneous equation corresponding 
to Eq. (|2^) has solutions and hence Eq. (|26| ) cannot be solved directly. We again intro- 
duce regularisation which enables a straightforward numerical solution. Thus we solve the 
equation 

2A /"°° 

k^{p, k) = M{p, k) + — P dq M^{p, q)K^{q, k), (28) 

TT Jo 

where 

M^(p,g) = M(p,g) + -^^. (29) 

For a sufficiently large cutoff A, the solution K{p, k) = K^{p, k) satisfies Eq. (|26| ) with very 

good accuracy and has oscillating behaviour for large p. K^{jp,k) is also oscillating with 

respect to A. Just like the K^ of the original problem, the solutions K'^ have the property 

that the ratio 

k^^{p,k) - k^^{p,k) 

k^3^p,k)-k^^{p,k) 

does not depend on the momentum p. The amplitude structure is given accordingly by 

K^'ip, k) = Koip, k) + A(A)ir,(p, k) (30) 

where Kq{p, k) is a solution to Eq. (p^) and Kh{p, k) = K^^{p, k) — K^^{p, k) is a solution 



to the corresponding homogeneous equation. Everything is the same as before except that 
now one can actually solve the homogeneous equation analytically. We find that there 
are two linearly independent solutions, sin [soln(a2p)] and cos [soln(a2p)], even though only 
one solution enters with a A-dependent coefficient in the general numerical form given by 
Eq. (|30|) . It is easy to see that the other linearly independent solution contributes into 
Ko{p, k) and is responsible for its oscillating behaviour. 

B. Extracting the physical ampUtude 

In the previous subsection we have shown that the inhomogeneous K matrix equation, 
Eq. ([T8|), can be solved numerically by introducing a sufficiently large cutoff A, and that the 
resulting solution is of the form given by Eq. (EHI). Writing Eq. (Bl) as 



K''ip,k)=K,,ip,k) + Ki^{p,k), (31) 

we are left with the task of extracting the physical K matrix Kph{p, k) from the numerical 
values for K^{p, k). Because Kph{p, k) vanishes for large p, the function K^{p^ k) in Eq. ( |3TD 
is that solution to the homogeneous equation corresponding to Eq. ([18|) which has the same 
asymptotic behaviour as K^{p, k); that is, K^{p^ k) satisfies the two equations 
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Kl^{p,k) = 


2A 


Jo 
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q' 


aq ivi {p, q 

q2 


-P 




oo 


{p,k). 







<(g,fc), (32) 

(33) 

Note that K^{p,k) cannot be obtained by solving the homogeneous equation correspond- 
ing to Eq. (pO|); this homogeneous equation has no non-trivial solutions since the inverse 
operator (1 — M^Gq)~^ exists. On the other hand, since Eq. ( |32D has solutions for every 
energy, its numerical solution presents a significant technical problem. For example, if one 
simply discretises this homogeneous equation and then tries to solve the resulting simulta- 
neous equations, one encounters severe numerical difficulties coming from the fact that the 
eigenvalues of the operator (1 — MGq) include not only zero, but also values infinitely close 
to zero. 

To solve Eq. (|3^) , and moreover, to obtain the solution with the asymptotic behaviour 
required by Eq. (|33|), we solve a sequence of inhomogeneous equations, like Eq. (|20|), but 
with progressively smaller inhomogeneous terms. That is, we solve the equations 

2A f°° a^ 

K^{p,k) = M,{p,k) + —P dq M\p,q)^^K^iq,k), (34) 

IT Jo q^ - k^ 

where i = 1,2,3,... and Mi{p,k) is taken to be a decreasing sequence of functions. In 
particular, we take 

Since Mi{p,k) becomes a vanishingly small function for increasing values of i, we expect 
K^{p, k), with A large enough, to approximate the solution of Eq. (|3^) for sufficiently large 
values of i. At the same time, since Mi{p,k) vanishes as p ^ oo, the solution K^{p,k) 
might be expected to have, up to some normalisation, the same asymptotic behaviour (as a 
function of p) as K^{p, k). That is in fact what we find empirically, so that 

Kt{p,k) ,^ a,Kt{p.k) (36) 

«— >oo 

where Cj is a constant. 

It is useful to illustrate the above procedure on the case of the homogeneous equation 
corresponding to Eq. ( pB]) for which an analytic solution is available. Thus we would first 
like to consider the equation 

2A /"°° 

kt{p. k) = M,{p,k) + —P dq M\p, q)K^{q, k), (37) 

TT Jo 

where Mj(p, k) is given by Eq. ( p5| ) and A has a large fixed value. By the above discussion, 
we expect the solution to Eq. (pT]) to have the form 

K^ip, k) = AKJ^ip, k) + i?f (p, k) (38) 
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that have the same 



FIG. 3. Solutions to the homogeneous equation corresponding to Eq. (| 
asymptotic behaviour as the solution to Eq. (|2^ ) with A 
solution KI^ to the homogeneous equation, the dot-dashed, dashed, and dotted lines correspond to 
the sequential approximations K^ with i = 4, 8, and 10, respectively. 



1.6 X 10^ Og . The solid line is the exact 



where K^{j), k) is that solution of the homogeneous equation corresponding to Eq. (p6D 
which has the same asymptotic behaviour as the solution K^{p, k) of Eq. (pSj), and Rf{p, k) 
is a remainder term that becomes progressively smaller with increasing i and that decreases 
rapidly as a function of p beyond a close neighbourhood of the origin. The constant Ai 
is expected to depend only on the choice of functions Mi{p,k). To see how these expec- 
tations are borne out in practice, we have numerically solved Eq. ( ^Tf ) for the sequence of 
inhomogeneous functions given by Eq. (|35D . We find that for different i and up to a normal- 
isation factor, the functions K^{p, k) all have identical oscillating behaviour for large p. We 
scale these functions so that their oscillating tails have the same amplitude, i.e., we rescale 
Eq. ( ^81) as 



K^{p,k) = kl^{p,k) + Rt{p,k) 



(39) 



with K^ = K^/Ai and Rf = R^/Ai, and then plot the results in Fig. ||. We see that in 
agreement with our expectations for increasing values of i, the remainder term Rf- becomes 
progressively smaller in such a way that two sequential functions K^{p,k) and K^_^^{p,k) 
coincide beyond some value of p that is progressively approaching zero. According to Eq. ( |39[ ) 
this tail, which is common for all functions K^, is just the tail of the solution to the 
homogeneous equation, K^. Choosing a different sequence of inhomogeneous terms, for 
example 



Mi{p,k) = a2 exp[-i{a2pY 



(40) 



we obtain a sequence of functions Kj^' that approaches exactly the same solution Kj^ to the 
homogeneous equation. For this non-physical case we are able to solve the homogeneous 
equation analytically. As the general solution is a linear combination of the independent 
solutions sin[soln(a2p)] and cos[soln(a2j9)], we are able to write our sequence of solutions as 



K^{p, k) = Ai sin[so ln(a2p) + a] + R^{p, k). 



(41) 
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FIG. 4. Numerical construction of the function Kj^{p,k) defined by Eq. (13 2| ) and Eq. (| 
The solid line is the solution K {p,k) to the inhomogeneous K matrix equation with cutoff 



A 



1.6 X lO^Og and on-shell momentum k 



0.82ao 



The dot-dashed, dashed, and dotted 



curves correspond to sequential approximations to KHp, k) with i = 4, 8, and 10, respectively. 



It is interesting to note that although the constant a depends on the cutoff A, our numerical 
findings indicate that it does not depend on the value of i. The function sin[soln(a2p) + a] 
(scaled appropriately) is also plotted in Fig. ^ As we can see, by using the sequence of func- 
tions approach, we can reproduce the particular exact solution to the homogeneous equation 
with better and better accuracy. Finally, by solving the inhomogeneous equation with cut- 
off, Eq. (pSD , we find that the asymptotic behaviour of the solution K^{p, k) is identical (up 
to a normalisation factor) with the asymptotic behaviour of all the solutions K^{p,k) of 
Eq. (pTl). It is thus K^{p, k) to which the sequence of functions K^{p, k) converges. 

We now apply the above procedure to our physical case of boson-dimeron scattering 

with a sufficiently large value of A in 

with the same value 

I5f). As in the simpler 



described by Eq. (|T8|). Thus we first solve Eq. ( 

order to obtain a particular solution K^{p, k). Then we solve Eq. ( 

of A using the sequence of inhomogeneous terms specified by Eq. 

case above, the asymptotic tails of K^{p,k) and all the solutions K^{p,k) are oscillatory 

and, up to a normalisation factor, identical. This can be seen explicitly in Fig. ^ where the 

Kf" for i = 4, 8, and 10 are plotted after having been scaled to have the same asymptotic 

tail as K^. From this it is seen how the sequence of functions K^ provides successively 

better approximations to the function K^ appearing in Eq. (|3TD . With i sufficiently large, 

we can subtract K^ from K^, in this way deducing the desired physical K matrix Kp^ with 

vanishing asymptotics. 

As evident from Fig. |^, lager and larger values of i are needed in order to approximate 
the value of K{^{p, k) for progressively smaller values of the momentum p. For our purposes 
it is sufficient to have i only large enough to give an accurate approximation to Kj^^p, k) at 
the on-mass-shell value of p, i.e., at p = k. By repeating the above calculations for a range 
of values of the on- mass-shell momentum, we obtain the physical K matrix Kph{k, k) as a 
smooth function of k. The result is drawn in Fig. ^. 

The above numerical procedure cannot be used directly to obtain the boson-dimeron 
scattering length because our sequence of functions method to determine Kj^^p, k) does not 
extend to the value of p = 0. However, taking into account the smoothness of the scattering 
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0.4 0.6 0.8 1.0 1.2 

k [l/aj 
FIG. 5. The on-mass-shell K matrix Kpfi{k, k) corresponding to the physical amplitude for 

boson-dimeron scattering. This numerical result is based on leading order EFT with two-body 

input only. 



amplitude, we can determine the scattering length by extrapolating the amplitude to zero 
momentum. There is not much sense in making that kind of extrapolation for the present 
scalar particle EFT model; however, we will be able to compare the predicted value for the 
scattering length with the experimental value when we complete calculations for the case of 
neutron-deuteron scattering in the doublet channel. 

IV. SUMMARY 



Applying leading order EFT to describe the scattering of a boson off a two-boson bound 
state results in an unusual inhomogeneous integral equation, Eq. (|^), originally derived in the 
1950's by Skorniakov and Ter-Martirosian |jll|, whose kernel MGq is such that the operator 
1 — MGq has no inverse. The same equation arises in the application of leading order EFT 
to neutron-deuteron scattering where its unusual properties have caused difficulties with its 
numerical solution |^. The problem encountered in Ref. ^ is that the scattering amplitude 
has no limit as the cutoff, used to solve the equation, is removed; in particular, it was found 
that as the cutoff is taken to infinity, the amplitude oscillates with a constant amplitude. 
It was argued that this undesirable behaviour of the amplitude should be eliminated by the 
introduction of a carefully chosen three-body force into the leading order Lagrangian ^ . 

In the present paper we have solved the problem of the ultraviolet cutoff-dependence 
of the scattering amplitude without the need to introduce a three-body force. We did this 
by taking into account the fact that, for any given energy, Eq. (^) has an infinite number 
of solutions amongst which only one corresponds to the physical amplitude for particle- 
bound state scattering. Rather than solving Eq. (P) directly, we chose to work with the 
corresponding K matrix equation, Eq. (plsl) . Although Eq. (|l^) still has an infinite number 
of solutions, it has the advantage of singling out only those amplitudes which are unitary. 



By introducing a cutoff A into Eq. (|T8D and numerically solving the resulting equation, 
Eq. (^), we have shown that every A that is sufficiently large, results in an amplitude that 
satisfies the original equation, Eq. (10). Thus, the sensitivity of the scattering amplitude to 
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the ultraviolet cutoff A is simply a reffection of the fact that one is in this way obtaining 
different solutions of Eq. (P) for every A. A simple numerical analysis of our solutions to 
Eq. (^) has allowed us to deduce the functional form of the A dependence of the solution - 
see Eq. (H). 

Amongst the infinite number of solutions K^{p, k) described by Eq. (p2D, the K matrix 
corresponding to the physical amplitude, Kph, is distinguished by the fact that it vanishes 
in the limit of infinite off-shell momentum p. This asymptotic behaviour can be contrasted 
with the oscillatory asymptotic behaviour of all other solutions. By introducing a sequence 
of equations like Eq. ( PO] ) but with progressively smaller inhomogeneous terms, we have 
managed to extract Kph numerically from our numerical solutions of Eq. (^0|) - see Fig. |^. 
In this way we have shown explicitly that EFT in leading order describes the scattering of 
a particle off a two-body bound state in terms of two-body forces only. 
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